General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



r 


1 





;X-.- 

' 

b V 

P v 5 * . 

V t ■1-. 

• - 

•; | 

■ 


; * 

-^r _>; ■ 

: * ‘ £ • 


• -- ' 

. -r • >• ■ 4 . 


VARIATION OF PARAMETERS FOR 
LONG PERIODIC TERMS IN PREDICTING 




LIFETIMES OF VENUS ORBITERS 

. . ■, - - r - ^ . > v ’ — - . /> ; ' 


V' v 

■* A ‘ ‘ 


- . «■ 

». 

• • « : . . ' > 

T ‘ : ■ “ . * / • 
7T - 'V > ' 


BERNARD KAUFMAN 


Wift ' 


.mm/. 

ife. 

rm* - ( 

g ' 

'» 


,. v •• ; f;f 

SEPTEMBER 1969 


GODDARD SPACE 

GREENBELT 

1 - / ■ * * 

“ir 


(THRU) 


(PAGES) 


(CATEGORY) 




X-531-69-395 

PREPRINT 


VARIATION OF PARAMETERS FOR LONG PERIODIC TERMS 
IN PREDICTING LIFETIMES OF VENUS ORBITERS 


Bernard Kaufman 


September, 1969 


Goddard Space Flight Center 
Mission and Trajectory Analysis Division 
Mission and Systems Analysis Branch 






^ 1 • 
* - 
V 4 



7 


PRECEDING PAGE BLANK NOT FILMED. 




VARIATION OF PARAMETERS FOR LONG PERIODIC TERMS 
IN PREDICTING LIFETIMES OF VENUS ORBITERS 





Bernard Kaufman 


f , ABSTRACT 

A very rapid variation of parameters method has been programmed to de- 
termine lifetimes of satellites in orbit about Venus. The perturbing function con- 
tains only the gravitational potential of the sun. The solar potential is expanded 
in a trigonometric series and then averaged over not only the period of the satel- 
lite but also over the period of Venus. The equation for variation of eccentricity 
uses only the singly averaged function so as to include medium period effects, 
however, for inclination and argument ofpericenter the doubly averaged function 
is used. 



It appears from the results of several sample cases that not only is this 
method considerably faster than a numerical integration solution of the three 
body problem (less than 15 seconds as opposed to one hour of machine time for 
an orbit study of 500 days) but is also highly accurate. Although this method 
contains only third body terms it should be very useful in first approximation 
studies. 
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VARIATION OF PARAMETERS FOR LONG PERIODIC TERMS 
IN PREDICTING LIFETIMES OF VENUS ORBITERS 


INTRODUCTION 

One of the determining factors in considering a Venus orbiter is the length 
of time the satellite will remain above the dense atmosphere of the planet. The 
orbit must be chosen so as to allow a sufficient lifetime to conduct tha required 
scientific experiments, but also to meet stringent quarantine conditions that 
have been imposed for Venus. 

In the present study, the atmosphere and oblateness of Venus is not con- 
sidered in the equations of motion. It is assumed that the gravitational attrac- 
tion of the sun is the only disturbing force acting on the satellite and that if the 
closest approach of the orbit is below a certain limit, atmospheric drag will 
cause the satellite to impact within a short time period. The radius of closest 
approach (periapsis r ) is defined as 


r p = a (! - e ) 


(1) 


where a is the semi major axis and e is the eccentricity of the orbit. It will be 
shown later that under long period perturbations due to the gravitational presence 
of a third body, the initial value of the semi major axis remains constant. There- 
fore by choosing a limiting value for r , equation (1) may be solved for a so- 
called critical value, e cr , of the eccentricity. A value of e larger than e cr then 
will be an impact orbit. 

The above outlined problem can of course be easily solved using a precision 
integration trajectory program to solve either Cowell's or Encke's equations of 
motion. However, past experience has shown that for small orbits, the computer 
time is prohibitive. This report investigates a variation of parameters solution 
that is intended to yield accurate but fast results. 

The Solar Potential Function: ^ 
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With reference to Figure A the solar potential IU may be defined as 



SUN 


K=P @ 


I- 111 

P~ r ' 3 


( 2 ) 


where 


and 


P - r - r 


_ _ f) 2 = r' 2 + r 2 - 2 r r' cos S 


or 


r = r 


" ( 1 + t 1 -^ cosS ) 


Then equation (2) becomes 


_ ^0 

/. r 2 2 r „ ^ 

1 + — COS S 

-1/2 

, r cos S 

r' 

_\ r' 2 r' > 

' r' J 


(3) 
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In this equation we may assume that r/r' « 1 and expand the square to include 
only terms of order 2. We then obtain 



In the usual form of the equations of motion the gradient of the above equation 
would yield the disturbing function. However, the gradient is taken with respect 
to the state of the sateilite and the first term of equation (4) contains only ele- 
ments of the sun. Therefore we have 



and we can redefine the solar potential as 


or 



( 5 ) 


Letting = n' 2 a' 3 where n' is the mean motion of the sun and a' 3 is the 
semi-major axis we can write equation (5) as 


R 


a 



[3 cos 2 S - 1 ] . 


( 6 ) 
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Now 


cos S = 


r r 


= r 


and 


?° = P cos d + Q sin 6 


( 7 ) 

( 8 ) 


where t) is the true anomaly and 


P = 


cos fl cos tu - sin fl sin w cos i 
sin fi cos a; + cos fl s in oj cos i 
sin i sincu 


Q = 


-cos flsin a' - sin flcos wcos i 
- sin flsin u + cos flcos wcos i 
s in i cos u> 


O) 


P is a unit vector in the satellite's orbital plane from the center of Venus to 
pericenter, Q is a unit vector in the orbit plane j)erpendicular to P in the direc- 
tion of satellite motion. The position of the sun in Venus' orbit plane may be 
written as 

( cos fl 
sin fl 
0 

Appendix A describes a venus centered coordinate system where the fundamental 
plane is Venus' orbital plane and also develops the transformation from the 
ecliptic plane to the new plane. 
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Equation (7) may now be written as 


denoting 


cos S = (P • r'°) cos 9 + (Q • r'°) sin 9 


P • r'° = a 


0 • r'° = 0 


cos S = a cos 9 + R sin 9 


where it should be pointed out that a and R are independent of S . Substituting 
equation (12) into equation (6) 


n * 2 /r \ 2 /£_\ 3 

2 U ) \ r ' ) 


[3 (a 2 cos 2 0 + 2a/3sin0cos0+/3 J sin 2 9) - 1] 


which after some manipulations may be written as 


R ® 2 


r — (a) ja 2 + |/3 2 - 1^ + i(a 2 - ( fl 2 )cos 2(9 + 3a/^sin 2(9j ( 13 ) 


Since we are mainly interested in long period terms it would be advantageous 
to eliminate any dependence of the above disturbing function on the mean 
anomaly i. This is done by means of the averaging process where the equation 
is integrated over the period of the satellite. The details of this process may 
be found in reference 1 and only the results are presented here. Equation (13) 
contains three terms which are dependent on the anomaly and the averages of 
these terms are: 
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Substituting the above equations into equation (13) we obtain 



a result which is considerably easier to use than is equation (13). 
Variation of Parameters 


( 14 ) 


( 15 ) 


(16) 


(17) 


The development of equations for the method of variation of parameters may 
be found in any good textbook on celestial mechanics and therefore will only be 
listed here. For completeness all six equations will be listed although it will be 
shown that not all of them are needed. 


d a 
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y \ — e 2 ° R Cr> 

dT " 
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dT~ 
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dR ® 
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Variation in Eccentricity: Long and Medium Period Terms 

The first of the above equations yields an extremely useful result. Since we 
have elliminated any dependence of R upon i we have 


= 0 
dt 


(19) 


which means that there are no long term or secular variations in the semi major 
axis. This means we need only to compvle the variations in eccentricity in order 
to obtain a time history of periapsis. We therefore need to look at the second 
of equations (26) and note that again 


3R„ 


34 


= 0 


and therefore 


de _ /l - e 2 
nea 2 


( 20 ) 


With reference to equation (17) we may write 


™ e Ba 


3R 


s 3/8 

3a 3 a) 3/3 3a> 


( 21 ) 


I 

» 


recalling the definitions of a, /3 , P and Q we have 


3q 3P _ 
3 oj 3 a> 


■i 

3 
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v-- 


and from equation (9) we have 


Therefore 


Similarly v/e find 


^=Q 


Q • r'° =/3 

d oj 


Then equation (21) becomes 


3 6) 3 a 3/3 


The derivatives are easily found to be 


?R, 2 r »2 /„* \3 

^ } (3 a) (1 + 4 e 2 ) 

3a 2 \ T . I 

3R„ „2 n <2 / < \ j 

-r# = I — } (3/3) (1 -e 2 ) 

a/3 2 Sr . y 


Substituting (25) into (24) and then into (2C) we obtain, after some algebraic 
manipulations: 


Je _ /r~2 15 


= - /\ - e 1 - - 
dt 2 
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where 

n' 

7 = — 

n 


The above equation represents the variation in eccentricity due to the 
gravitational presence of the Sun. In the case of Venus, we have a near circular 
orbit for the sun and therefore we may simplify equation (26) further by setting 





Equation (26) also indicates that in order to completely solve the equation 
we also need a time variation of fi, co and i which appear in a and /3. Before 
developing these equations we will make a further simplification of R 0 which 
will reduce the system to one of only 3 coupled equations. 


Medium Period Terms of R . 

Lorell and Anderson (Reference 2) have made a further approximation by 
averaging over the period of the third body about the central body. In our case 
this would mean eliminating dependence upon fi 0 defined earlier as the longitude 
of the sun measured in Venus' orbital plane. From equation (17) and the defini- 
tions of a and /3 we need only consider the two terms 


and 



a 2 - /3 2 


(27) 

(28) 


[ 



we consider first (27) 

a 2 + /3 2 = cos 2 fi cos 2 fi 0 + sin 2 fi cos 2 i cos 2 fi 0 

+ sin 2 fisin 2 fi 0 + cos 2 ficos 2 i s in 2 fi 0 
+ 2 sin fl cos fi sin fi 0 cos fi 0 
- 2 s in fi cos fi cos 2 i s it fi Q cos fi Q 
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= 1 -sin 2 H 0 cos 2Q sin 2 i + sin 2 i 


1 s i n 2f2sin 2 Q_. 

2 ® 


-sin 



s in 2 i 

2 


[ 1 - cos 2 fl cos 2 fi Q ] 


+ 


s i n 2 i 

2 


Lsin 2ft s in 2 ft 0 


1 - 


sin‘ 


l s in* 
— + — 


— cos 2 (ft - ft) 


Therefore 


3 

2 




3 . 2 . 

T- Sin 1 

4 


3 • 2 

+ — s in' 

4 


i cos 2 (ft 0 - fl) 


(29) 


In a similar manner we find that 


a 2 - <3 2 = — sin 2 i cos 2a; + cos 4 ( -i- ] cos 2 (f2 0 - a) - fi) 


a) 


+ S 


- (I) 


cos 2 (Q + c-j - H) 


(30) 


when we carry out the averaging process we find that all terms containing f2 0 
may be dropped. We may therefore redefine equation (17) as 


R© = a 2 n ' 2 


/ 1 3 A (. 3e 2 \ 15 e 2 . 2 • 0 . 

_\4 8 ) \ 2 ) 16 K J 


( 31 ) 


where we have let 1- = i 
r ’ 
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Variation in Eccentricity: Long Periodic Terms Only. 

With the new definition of R 0 above we return now to the second of equations 
(18) and calculate again de /dt. Again we need only calculate BR @ /Boj 


1 c 

® 2 #2 2 • 2 * • o 

— — a n e sin* i sin 2 oj 

8 


de (1 - e 2 ) 1 2 f 15 2 <2 2 • 2 • 0 1 

dt pnfl 2 L 8 J 


de 15 n' 2 2 xi /2 ■ 2 • • o 

— — = — e (1 - e , ')* / '' sin'' i sin 2u> 

dt 8 n ' 


This equation is considerably simpler than the previous equation (26) and it 
is noted that the dependence upon fi has been eliminated. However, the depend- 
ence upon i and o> is still present and we must now develop the equations for 
these terms. 


Variation in Inclination and Argument of Periapsis 

Using the new definition of R 0 above we return now to the third of equations 


1 P R 0 . 9 R o" 

TT7 2 sin i j 


We observe immediately that 


a result that was not true before the elimination of medium period terms. Also 
from equation (32) 
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15 2 ' 2 2 2 • o 

= a n e sin z i sin 2 go 

du) 8 


substituting this we obtain 


or 


di 

dt 


15 

8 n 



s in i cos i s in 


2 co 


di 

dt 


15 

16 n 


- sin 2i sin 2 co 

/l - e 2 


(34) 




The equation for the variation in argument ofperiapsis requires dR 0 /Be and 

/di. 



15 , . 2 

+ -g-( sin 


i cos 2w) 



(s in i cos i cos 2co) 


Substituting these into the fourth of equations (18) 


d a) 

d t 


<2 

H_(l - e 2 ) 1/2 



9 . 2 . 15 . . 2 

— sin' 4 l + — (s in' 4 
8 8 


i cos 2w 


) 


n' 2 (1 - e 2 ) 1 2 

n (1-e 2 ) 




3 2 . 

— COS' 4 1 

4 


(. 3 e 2 \ 

15 2 2 • o 

1 + 

+ — c 4 cos' 4 i cos 2 oj 

' 2 y 

8 
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• -f: ■ , • - .?;>• ;-' -'f' -‘ • ‘ - ’ *•>„■ , 


- — (1 - e 2 ) 1/2 J ^-- |-sin 2 i (3 - 5 cos 2o>) 


+ i cos2 1 + 1 e2 cos2 1 _ H e 2 COS ' 1 cos 2a) 
4 (1 - e 2 ) 8 (1 - e 2 ) 8 (1 - e 2 ) 


_ n ' 2 (1 - e 2 ) 1/2 J 3 3 . 


sin 2 i (5 sin 2 co - 1) 

4 4 


3 cos 2 i 3 e 2 cos 2 i 
+ 4 (1 - e 2 ) + 4 (1 - e 2 ) 


(5 s in 2 u> - 1) 


— - — (1 - e 2 ) 1/2 <1-5 sin 2 i sin 2 oj + sin 2 i 
4 n 


1 - 5 e 2 sin 2 a>-e 2 e 2 sin 2 i - sin 2 i - 5 e 2 sin 2 i sin 2 cv 
(1-e 2 ) (1 -e 2 ) 


3n' 2 2 \i/2 fc -2 fe 2 - e 2 s in 2 i • 2 -1 

4 " \ L 1 -e 2 J 


+ 1 + s in 2 i 


1-e 2 e 2 sin 2 i - sin 2 i 
+ + 

1-e 2 1-e 2 


and finally 


*2L = 1^1(1 - e 2 ) 1/2 Jl + 

dt 2 n 1 


Isin^ (ei - Sin2i) \ 

2 (1 - *’) J 


Equations (33), (34) and (35) represent a system of three equations in e , i 
and a; only and it is now possible to solve for a time history of the eccentricity 
without having to solve for % or fl. Actually one may use equation (26) instead 
of equation (33) by merely holding ft constant atn =fi Q . in the computer program, 
equation (26) was utilized instead of (33). The results will not differ significantly. 
It is to be noted that the above three equations are identical to those given by 
Lorrel and Anderson in reference 2. 
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Results 


Figures 1 through 5 show comparisons between the above outlined method 
and a precision n body Encke integration program. As can be seen from these 
graphs there are no significant differences between the results obtained. How- 
ever, it must be pointed out that the n body program took approximately one hour 
to integrate the orbit for 500 days and the variation of parameters method took 
less than 15 seconds. 

Figures 1 through 3 are of similar size orbits with only the inclinations being 
changed. The apparent secular variation in eccentricity appears in all three 
graphs, however for an inclination of 5°, this variation is very small. Figures 
4 and 5 show trajectories that impact the planet after 500 and 217 days respec- 
tively. Impact is defined as that point where the Venusian atmosphere will 
drastically alter the orbit such that virtual impact will shortly occur. For this 
study an altitude of 200 km was used for impact. 

Figure 6 through 9 '•.re for orbits with periapses of 6575 km and apoapses of 
46025 km. Here the parameters that were varied are inclination and argument of 
pericenter. All four of these orbits impact in less than one terrestial year. 

Figures 10 through 14 represent an orbit with periapsis radius of 7550 km 
and apoapsis radius of 46050 km. Again the parameters that were varied are 
inclination and argument of pericenter. The first three of these orbits have not 
yet impacted after 800 days but all indications are that they will. The last two 
will impact within a very short time after 800 days. 
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APPENDIX A 


TRANSFORMATION FROM ECLIPTIC TO VENUS' ORBITAL PLANE 



Figure A-l-Venus Coordinate System 


Define: 

i ? = inclination of Venus' orbital plane to ecliptic 

= right ascension of ascending node of Venus' orbit upon the ecliptic 


29 


fi t = right ascension of ascending node of satellite orbit upon ecliptic 

a> t = argument of periapsis of satellite orbit referenced to ecliptic 

N ( = ascending node of satellite referenced to ecliptic 

fl = right ascension of ascending node of satellite orbit upon orbital 
plane of Venus 

oo = argument of periapsis of satellite orbit referenced to orbital plane 
of Venus 

N ? = ascending node of satellite referenced to orbital plane of Venus. 

Let be the Venus centered position of the Sun in mean ecliptic and equii.ox 
of Date Jan. 1, 1972 at 0*’ r . 



/cos 

D' 

cos 


/ 

© 


11 
© o 

II 

I s in 


cos 


V 

s in 

i>' 



(Al) 


where at the above date 

D' 0 = 175°2756: - 76°427751766 

<*>' 0 = 3°3542; i 9 ^ 3°394354251 
Let 

T = *..1 (ij) *..i < n 9> (“) 


denote a rotation first about the z axis of the ecliptic through the angle £2 9 fol- 
lowed by a rotation about the new x axis of the ecliptic through the angle i ? . T 
is then the transformation matrix from the ecl'ptic plane to the orbital plane of 
Venus with the x axis in the direction of the ascending node of Venus’ orbit upon 
the ecliptic and the Y axis 90° in the direction of planetary motion. 
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K • 


Then r° 0 , the position of the Sun in Venus' orbital plane at the above date, 


is 


r 

K 

§ 

s 


also 


(1972) = tan" 1 


^ cos 


To obtain at any time t we must have then 

n; t =n* @ (1972) + 1.6 At 


(A3) 


(A4) 


m 

R 

€ 




S. 


where At is days from Jan. 1, 0 hr • 1972. Equation (A4) is needed if the form of 
the disturbing function in equation (17) is used. 

The orbital elements of the satellite referenced to Venus' orbital plane may 
now be found as follows: 

The position r° of the satellite referenced to the ecliptic is 

r° = P* cos 8 + Q' sin 8 


where 


P' 


cos fl , cos a) , - sin A , sin oj , cos i , 

S€ S€ s€ S€ 

s in fi , cosoi , + cos fl , s in a> cos i , 

S« S( s€ Sf Sf 


Q' = 


sin k s in oj 
sf se 


cos A , sinco,-sinfl , cos co , cos i , 

S€ s€ S € S€ s€ 


sin fl , siniu , +cos fl , cos w. , cos i. 

S € S € S € S€ S € 


sin l ^ cos co _ 

S € $€ 


(A5) 
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and 6 is the true anomaly of the satellite. Then f the position of the satellite 
referenced to Venus' orbital plane is 

r 0 = T r*° - (TP') cos 8 + (TQ') sin 8 
Letting TP' = P; TQ' = Q 


r° = P cos 8 + Q sin 9 


(A6) 


where P and Q are defined as above in (A5) with the appropriate change in nota- 
tion to signify that the reference plane is now the orbital plane of Venus. Using 
the notation 


P = 



Q = 



and 


we have 



(A7) 


and 


i s$ = cos 1 (Pj Q, - P 2 Qj) = cos 1 (R 3 ) 


fi - Tan 1 
s 9 


/ P 2 Q, - P 3 Q \ 
\ P i Q 3 - p 3 Qij 


= Tan" 1 



(A 8) 


(A 9) 
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Finally we write the full expansion for 


#*r, 

> 


T = 


cos sin fly 

sinfi^cos i,j cosfl^cosi^ 
sin n 9 sin i 9 - cos sin i 9 


0 

sin i^ 
cos i^ 


(A10) 
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